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Abstract. We present radio and submillimetre observations of the 04I(n)f star £ Pup, and discuss structure in the outer region 
of its wind (~ 10 - 100 R,). The properties of bremsstrahlung, the dominant emission process at these wavelengths, make it 
sensitive to structure and allow us to study how the amount of structure changes in the wind by comparing the fluxes at different 
wavelengths. Possible forms of structure at these distances include Corotating Interaction Regions (CIRs), stochastic clumping, 
a disk or a polar enhancement. As the CIRs are azimuthally asymmetric, they should result in variability at submillimetre or 
radio wavelengths. To look for this variability, we acquired 3.6 and 6 cm observations with the Australia Telescope Compact 
Array (ATCA), covering about two rotational periods of the star. We supplemented these with archive observations from the 
NRAO Very Large Array (VLA), which cover a much longer time scale. We did not find variability at more than the +20 % 
level. The long integration time does allow an accurate determination of the fluxes at 3.6 and 6 cm. Converting these fluxes into 
a mass loss rate, we find M = 3.5 x 10~ 6 Mq/jx. This value confirms the significant discrepancy with the mass loss rate derived 
from the Ha profile, making f Pup an exception to the usually good agreement between the Ha and radio mass loss rates. To 
study the run of structure as a function of distance, we supplemented the ATCA data by observing £ Pup at 850 fjm with the 
James Clerk Maxwell Telescope (JCMT) and at 20 cm with the VLA. A smooth wind model shows that the millimetre fluxes are 
too high compared to the radio fluxes. While recombination of helium in the outer wind cannot be discounted as an explanation, 
the wealth of evidence for structure strongly suggests this as the explanation for the discrepancy. Model calculations show that 
the structure needs to be present in the inner ~ 70 R, of the wind, but that it decays significantly, or maybe even disappears, 
beyond that radius. 
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1. Introduction 

Zeta Puppis (HD 66811; 04I(n)f) is one of the most-studied 
O-type stars. Its stellar wind is driven by radiation pressure, 
whereby the momentum of the stellar photons is transferred 
to the wind material via line opacity. With other early-type 
stars it shares the property that its wind is structured. Evidence 
for this structure for O-type stars in general is seen in the 
Discrete Absorption Components (DACs) and black troughs 
in the ultraviolet resonance lines (e.g. Prinja et al. 1990), the 
presence of X-ray emission (Sciortino et al. 1990 1 and the ex- 
cess flux at infrared and millimetre wavelengths (Runacres & 
Blomme 1996). Specifically for £ Pup, the significant discrep- 
ancy between the mass loss rate derived from the Ha spec- 
tral line and that from the radio continuum flux (Petrenz & 
Puls[T996 1 also points to structure in the wind. 

Various types of structure have been proposed to explain 
the observations. Inhomogeneities on the stellar surface (e.g., 
due to non-radial pulsations), or magnetic fields, result in 
somewhat different radiative forces, thereby creating fast and 
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slow streams of gas. As these streams also rotate, they col- 
lide, creating Corotating Interaction Regions (CIRs). These are 
large-scale, spiral-shaped structures in the wind that corotate 
with the stellar surface. In some cases, the CIRs can explain the 
most interesting property of DACs: the recurrence time scale is 
in agreement with the estimated rotation period (Mullan 1986 
Cranmer & Owocki [T996l . 

Another type of structure that can exist in the wind is of a 
smaller-scale, stochastic type. It is caused by the inherent in- 
stability of the radiative driving mechanism (Owocki 2000). 
Stochastic structure has been used to explain the X-ray fluxes 
(Lucv ll9^2al Hillier et al. lT9951 Feldmeier et al. 113971 . the 
black troughs in the UV resonance lines (Lucy 1982b ) and the 
excess millimetre flux (Blomme et al. 2002). 

Finally, it should be noted that f Pup is a rapid rotator (at 
43 % of critical velocity). This suggests the possibility that 
f Pup has a disk, or at least some enhancement of its density 
near the equatorial plane. Harries & Howarth (1996) derived 
a density contrast of at least 1.3 from their linear spectropo- 
larimetry of Ha. Such equatorial density enhancement has been 
used (Petrenz & Puis 1996) to explain the discrepancy that 
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smooth and spherically symmetric models give for the mass 
loss rates derived from Ha and radio observations. 

While observational evidence for the existence of disks in 
certain stars (e.g. Be stars) is strong, there are theoretical dif- 
ficulties in understanding how these disks form (Owocki et 
al. Q998). One of the explanations proposed was the wind- 
compressed disk model (Bjorkman & Cassinelli 1993 1, but a 
more detailed study showed that the non-radial component of 
the line force and gravity darkening can combine to inhibit the 
formation of the disk. Interestingly, in some circumstances, the 
detailed models show an enhancement of the density towards 
the poles rather than the equator (Owocki et al. 119961 [1998 
Petrenz & Puls l2000l 

Observations at submillimetre and radio wavelengths are 
relevant to structure. This is because bremsstrahlung, the domi- 
nant emission process at these wavelengths, has two interesting 
properties. First, its opacity is proportional to the wavelength 
squared. Therefore, the observed emission originates above a 
characteristic radius that increases with wavelength (typical 
values are ~ 5 R m at 1 mm and ~ 100 R* at 20 cm). Secondly, 
free-free emission also depends on the density squared. This 
makes it a good indicator of structure: a structured wind will 
have more radio emission than a smooth one. By looking for 
variability at a certain wavelength, one can hope to distin- 
guish between structure that has azimuthal symmetry (stochas- 
tic, disk or polar outflow) and structure that lacks azimuthal 
symmetry (CIRs). By looking at increasing wavelengths, one 
can see how the "amount" of structure changes in the wind. 

We have already applied these techniques to the submil- 
limetre and radio observations of the BOIa star e Ori (Blomme 
et al. I2002i . For this star, we did not detect variations in the 
radio emission (at the 25 % level), but we found that the mil- 
limetre flux is substantially higher than a smooth wind model 
predicts. This discrepancy was interpreted with a model for 
stochastic structure, showing that considerable structure must 
persist up to at least ~ 10 R t in the wind of e Ori. The present 
paper extends the study of structure in the outer wind to a lu- 
minous early O-type star. 

To look for variability at radio wavelengths, we acquired 
3.6 and 6 cm observations with the Australia Telescope 
Compact Array 1 (ATCA), during a 12-day observing session, 
which covers about two rotational periods of the star. To see 
how structure changes as a function of distance, observations 
were collected at 850 /mi with the James Clerk Maxwell 
Telescope 2 (JCMT), and at 20 cm with the NRAO Very Large 
Array 3 (VLA). We supplemented this material with data from 
the VLA archive. A log of all data used in this paper is given 
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Table 1. New observations and archive data of £ Pup. 



programme 


date 


A 


name 




(cm) 


New observations 


ATCA 






C824 


1999-09-15^27 


3.6+6 


VLA 






AB1017 


2002-03-27 


20 


JCMT 






M00BU20 


2000-10-12 


0.0850 


Archive data 


VLA 






FLOR 


1978-07-23 


6 


FLOR 


1978-10-13 


6 


BIEG 


1978-11-05 


6 


NEWE 


1979-02-09 


2 


FLOR 


1979-02-16 


20 


CHUR 


1979-07-12 


6 


FLOR 


1981-10-18 


1.3+2+6+20 


AA28 


1984-03-07 


2+6 


AB327 


1985-01-29 


2+6 


AH365 


1989-05-13 


3.6 


AC308 


1995-01-17 


20 



in Tabled Upper limits from other observations (not discussed 
here) are listed in Wendker ( 1995 1. 

In Sect. |2 we present the stellar parameters of £ Pup. The 
ATCA, VLA and JCMT observations are discussed in Sects. [3] 
|4] and [5] respectively. The interpretation of the observational 
material is discussed in Sect. [6] and conclusions are drawn in 
Sect.0 

2. The parameters of £ Pup 

Table |3 lists the £ Pup parameters we use in this paper. The 
stellar parameters and wind parameters are taken from the uni- 
fied NLTE model of Puis et al. (1996 1. These values are in 
close agreement with models by Bohannan et al. {I990l and 
Pauldrach et al. ( 1994i. The traditionally assumed distance of 
450 pc (Kudritzki et al. 1983 1 falls well within the Hipparcos 
error bar, so we will use d = 450 pc throughout this paper. The 
radius of 19 Rq at a distance of 450 pc corresponds to an an- 
gular diameter of 0.39 milli-arcsec (mas), in good agreement 
with interferometric measurements (#ld = 0.42 ± 0.03 mas; 
Hanbury Brown et al. U974i . 

£ Pup shows considerable evidence for structure in its 
stellar wind. Discrete Absorption Components (DACs) are 
seen to move through the ultraviolet resonance lines (Prinja 
et al. 1992 and references therein). In 1995, ( Pup was one 
of three targets observed continuously during 16 days by 
the International Ultraviolet Explorer (IUE Mega Campaign, 
Massa et al. 1995 1. The analysis of this high-quality dataset by 
Howarth et al. ( 1995 1 reveals a 19.2 h and a 5.2 day period. The 
19.2 h period is the recurrence time of the DACs. The 5.2 day 
period indicates modulation on a global scale, and is in good 
agreement with the rotation period estimated from v sin i, sug- 
gesting that material in the wind corotates with the star. One 



R. Blomme et al.: Radio and submillimetre observations of wind structure in ( Pup 



3 



Table 2. Stellar parameters of f Pup. 



parameter 


value 


reference 


RA (J2000) 


08 h 03 m 35 s :0467 


SIMBAD 


Dec (J2000) 


-40°00'11'.'332 


catalogue 


fi a cos S 


-31.7 ± 0.5 mas/yr 


G01 


LI s 


+ 17.6 + 0.6 mas/yr 


G01 


V magnitude 


2.25 


M87 


B- V 


-0.27 


M87 


^B-V 


0.044 


S77 


spectral type 


04I(n)f 


W72 




42000 K 


P96 




3.60 


P96 


I02 LI Lr7\ 


6.00 


P96 




19% 


P96 




52.5M Q 


P96 




0.12 


P96 


i' sin i 


220 km/s 


P96 


v«, 


2250 km/s 


P96 


/8 


1.15 


P96 


M 


5.9 x 10- 6 Mo/yr 


P96, from Ha 




450 pc 


K83 




429 + ™ pc 


Hipparcos 



References: 



G01 Gontcharov et al. <200"T1 

K83 Kudritzki et al. (19831 

M87 average from Mermilliod <1987t 

P96 Puis et al. fl99"6l 

S77 Snow et al. (19771 

W72 Walborn (19721 

possible cause for this corotating material is the presence of 
a (weak) magnetic field. Attempts to detect the magnetic field 
for ( Pup have so far yielded a null result with an uncertainty 
of 100-200 G (Barker et aL fTWil Chesneau & Moffat|l002 1. 

( Pup also shows variability in the Ha and He n 4686 
line profiles (Moffat & Michaud 1981 and references therein, 
Hendry & Bahng [TWA Reid & Howarth [19961 Eversberg et 
al. [H)98). The presence of variability in these wind-formed 
lines is consistent with large-scale structures corotating through 
the wind. Eversberg et al. however interpret their He n 4686 
observations in terms of smaller-scale structures ("clumps") 
moving out in the wind. Fullerton et al. ( 1996 > and Reid & 
Howarth ( 1996 1 detected variability in a number of optical lines 
(He i, He n, N iv, C iv). The variations have a period of 8.54 h. 
Both papers attempt to interpret the variations in terms of non- 
radial pulsations, but they also raise concerns that these lines 
might be wind contaminated. The photospheric perturbations 
due to non-radial pulsations are another possible way of creat- 
ing CIRs. 

Optical continuum fluxes also show variability: the 5.2 
day rotation period was detected in optical photometry 
(Balona 1992 1. From Hipparcos photometry, Marchenko et 
al. ([T998 1 derive a period which is half the rotation period. 

Other indicators of structure are the continuum flux ex- 
cess at infrared and millimetre wavelengths (Runacres & 
Blomme 1996) and the presence of X-ray emission (Long & 
White fT9"S0l . Models by Hillier et al. ( fT99"3l and Feldmeier 



et al. J1997I show that structure can explain the X-ray emis- 
sion. Earlier claims that the observed X-ray emission is vari- 
able (Collura et al. 1989 1 were later shown to be incorrect 
(Berghofer & Schmitt 1994). Later, Berghofer et al. (1994) 
found variability in the ROSAT X-ray data (on roughly the 
same time scale as the Ha variability), but this was not con- 
firmed by the ASCA data (Oskinova et al. 1200 lb . 

3. ATCA observations 

3.1. Data 

We used the Australia Telescope Compact Array (ATCA) for a 
long observing session on £ Pup in September 1999. The ~ 27 h 
on-target integration time was divided into 8 runs, spread over 
12 days, thus providing coverage of the -5.2 day rotation 
period. A detailed observing log is presented in Table [3] At 
the time of the observation, ATCA was in configuration 6A, 
with the longest baseline at 5939 m and the shortest one at 
337 m. The continuum observations were done simultaneously 
at 3.6 cm (X-band; 8.688 GHz) and 6 cm (C-band; 4.848 GHz). 
The observing bandwidth was 128 MHz. A single observing 
run consists of repetitively observing ( Pup for 15 min and 
then the phase calibrator PKS B0823-500 for 5 min. The flux 
calibrator PKS B 1934-638 was observed when possible at the 
beginning or end of the run. The flux calibrator could not be 
observed during the fourth (SEP20) and seventh observing run 
(SEP25) due to scheduling or technical problems. 

3.2. Reduction 

As the level of variability we expect is small, the reduction 
needs to be done as carefully as possible. For this reason, we 
also discuss the reduction in detail. The data were reduced in 
Miriad following the user guide (Sault & Killeen 1999). The 
data were read into Miriad and corrected for self-interference 
of the array, as well as for the phase difference between the X 
and Y channels. Next, the programme blflag was used to flag 
out bad datapoints interactively. The calibrators were then used 
to determine the antenna gains as a function of time (using mf- 
cal). 

The flux assigned to the flux calibrator PKS B1934-638 is 
5.85 Jy (with a 2 % error) at 6 cm. The flux value determined 
with uvflux of the phase calibrator PKS B0823-500 is 3.0 Jy. 
At 3.6 cm the flux calibrator is 2.88 Jy (with a 2 % error) and 
the phase calibrator is 1.53 Jy. In those runs where the flux 
calibrator had not been observed, we used the phase calibrator 
with the above mentioned flux values to calibrate the data. We 
determined the bandpass function from the flux calibrator. 

We interpolated the instrumental gains and applied them 
to the £ Pup observation. The task invert was then used to pro- 
duce an image from the visibility datasets by Fourier transform. 
In the inversion we used multi-frequency synthesis (MFS), 
which compensates for the spectral index of the source across 
the bandwidth. Because of the sparse UV-coverage per obser- 
vation, we used robust uniform weighting (Briggs 1995 ) to im- 
prove the RMS in the map. To de-convolve the image we used 
mfclean. The resulting clean components are then convolved 
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Table 3. ATCA data for f Pup. Start and end times of each observing run are given, as well as the total duration of the run and the 
integration time on f Pup. The measured fluxes at 3.6 and 6 cm are listed, with their error bars and the size of the beams (major 
x minor axis and position angle). The data obtained on SEP23-24 were split in two sets in order to have the same time span as in 
the other observing runs. 



start date-time 


end date-time 


total 


(Pup 


F 3 . 6 


cm 


r 6 cm 


beam 3 6 cm 


OA 

rA 


beam 6 cm 


1QQQ- 


1 QQQ- 
vyyy- 


time (h) time (h) 


(mJy) 








\ x ; 


oEr lj-Zj.4 / . JJ 


OJir 10-Uj. jZ. Lj 




O A A 

Z.44 


1 no _i_ 
l.Vo ± 


U. Iz 


1 ac\ -4- n 1 1 
1 .^u ± u. 1 1 


a nivn aq 

4 .UOXU.Oo 


jj.O 


7 AT, v 1 73 


ohr lo-z3:U /: 1 j 




3.52 


2.56 


Z.1U ± 


U. IZ 


1 ac\ 1 n 1 1 
1.4V ± U. 1 1 


3./8XU.O/ 


on a 
LyA 


A n/ivy 1 on 
0.9oXi.2U 


0r.rl7-zz.u7.UJ 


oCrZU-Ul .J7.JJ 


3.44 


2.50 


1.98 ± 


0.13 


1 .JH ± U, 1 1 


A 31 vO 77 


1 8 1 

1 O.J 


fi r»7 v 1 1 n 

O.U / A1.1U 


ccpon 1 
otirzyj- 1 o.UO.Zj 


oJ_,r ZU-ZZ. Ij. jj 


3.79 


2.92 


2.40 + 


0.12 


1 A1 -4- n 1 1 
1 .01 ± u. 1 1 


q ^Q v n A3 




A A3 v 1 1 A 


SEP22-20:08:55 


SEP23-02: 16:45 


5.26 


3.93 


2.22 ± 


0.13 


1.62 ±0.10 


2.65x0.62 


9.5 


4.82x1.14 


SEP23-15:07:55 


SEP24-01:38:55 


9.39 


6.86 


2.41 ± 


0.13 


1.67 ±0.11 


2.41x0.79 


-45.1 


4.35x1.44 










2.40 ± 


0.13 


1.66 ±0.10 


3.22x0.60 


9.2 


5.92x1.09 


SEP25-20:05:25 


SEP26-00:00:25 


3.54 


2.74 


2.55 ± 


0.13 


1.60 ±0.11 


4.11x0.57 


2.6 


7.74x1.04 


SEP27-16:04:05 


SEP27-19:59:35 


3.60 


2.62 


2.39 ± 


0.14 


1.53 ±0.11 


3.30x0.73 


-40.3 


6.02x1.32 


TOTAL 




35.84 


26.57 




















average flux 


2.26 




1.55 














median flux 


2.39 




1.60 












Complete data set 


2.38 ± 


0.09 


1.64 ±0.07 


1.65x0.73 


2.6 


3.01x1.33 



ATCA 1999-SEP 6 cm 

T I I I I I I I I I I I I I I I I I [■ 



S6* 

-39 C 5B' - '57 - 

"S3 

o 



-40°02' - ' S5 



j 1 1 I 1 1 1 1 I 1 1 , 1 I 1 1 1 1 I 

08 h 03 m 40 s 30 s 20 s 10 s 

RA (J2000) 

Fig. 1. f Pup and 6 other sources detected at 6 cm on the 
ATCA images. This map covers a large part of the primary 
beam. The synthesized beam is shown in the lower left cor- 
ner. Identifications of the sources are listed in Table 0] Sources 
S6 and S7 are not seen on the 3.6 cm image, because of the 
smaller primary beam. 



with a Gaussian and added to the residual image. Besides 
£ Pup, 6 other objects were detected in the primary beam at 
6 cm and 4 at 3.6 cm (see Fig.[2and Table0Ji. 

The cleaning was stopped when the absolute maximum in 
the residual map (i.e. the map from which the clean compo- 
nents have been subtracted) reached 0.30 mJy. This is about 3 
to 4 times the theoretical noise of each map, where the theoret- 
ical noise is calculated taking only the system temperature of 
the front-end receiver into account, not the calibration errors, 
side-lobes or any other instrumental effects. Experience shows 



Table 4. Position of other sources on the combined ATCA 
image. Formal error bars are better than 0TO1 in right as- 
cension and 071 in declination. Identifications given refer to 
Jones ( 1985 1 and Condon et al. ( 1998 1. No identification was 
found for S4 or S7 in the SIMBAD or HEASARC catalogues. 



no. 


RA (2000) 


DEC (J2000) 


Identification 




S2 


08 03 33.11 


-40 00 03.1 


EQ 0801-398 




S3 


08 03 28.65 


-39 58 20.4 


NVSS J080327- 


•395828 


S4 


08 03 24.38 


-39 59 49.2 






S5 


08 03 21.54 


-40 01 56.9 


NVSS J080321- 


•400153 


S6 


08 03 10.87 


-39 57 29.8 


NVSS J080310- 


•395729 


S7 


08 03 44.03 


-39 57 56.5 







that at this cutoff the RMS in the cleaned map and in the resid- 
ual map are both about equal to the theoretical RMS. At the 
end, to correct for primary beam attenuation, the task linmos 
was used. 

All images were reduced in exactly the same way with ex- 
actly the same parameters. 

3.3. Fluxes and error bars 

The flux of ( Pup was determined by fitting an elliptical 
Gaussian, with values for the major and minor axis and po- 
sition angle kept fixed to the beam values. Table [5]lists the flux 
values at 3.6 and 6 cm, and the respective beam sizes with their 
position angle. Because the observations at 3.6 and 6 cm were 
done simultaneously, the position angles are the same at both 
wavelengths. The data obtained on SEP23-24 were split in two 
sets in order to have the same time span as in the other observ- 
ing runs. The 6 cm fluxes of the other sources in the image are 
less than 1.2 mJy, except S2 and S6 which are about 4 mJy 
(uncorrected for the primary beam effect): these values are low 
enough that we do not have to worry about the effect of their 
sidelobes. The fitting procedure (imfit) gives an RMS error on 
the flux measurement. To take into account the calibration er- 
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ATCA 1999-SEP 3.6 cm 



date (1999-SEP-) 



-40°00'15" — C> 




0.80 0-90 1.00 1.10 1.20 1.30 
f Pup normalized flux 




0.80 0.90 1.00 1.10 1.20 1.30 
f Pup normalized flux 




0.70 0.80 0.90 1.00 1.10 1.20 1.30 
£ Pup normalized flux 




0.4 0.6 0.8 1.0 1.2 
C Pup normalized flux 



Fig. 2. The top panel shows the £ Pup fluxes of each ATCA ob- 
servation at 3.6 cm (O) and 6 cm (□) as a function of time. 
The flux derived from the single dataset combining all the 
ATCA observations is given by the dotted line. The lower 4 
panels show the correlation of the ( Pup fluxes with those of 
S2=EQ 0801-398 and S3=NVSS J080327-395828. In the cor- 
relation plots the fluxes were normalized by their average. 



08 03 35.5 35.0 s 34.5" 

-0.1,0.1,0.14,0.2,0.28,0.4,0.57,0.8,1.13,1.6,2.26 mjy 



ATCA 1999-SEP 6 cm 




-40°00'15" — 



08"03"'35.5° 35.0° 34.5° 

RA (J2000) 

-0.1,0.1,0.14,0.2,0.28,0.4,0.57,0.8,1.13 mjy 



Fig. 3. ATCA observation (all data combined) of f Pup at 3.6 
and 6 cm. The cross indicates the optical ICRS 2000.0 position 
(from SIMBAD), corrected for proper motion (see Table |2j. 
The contour levels follow a logarithmic scale: their values are 
listed at the bottom of each figure. The negative contour is 
given by the dashed line. The first positive contour is at about 
three times the RMS noise in the map. The beam is shown in 
the upper left corner. 



ror, we added 2 % of the flux values to this RMS, and thus 
arrived at the final error bars (listed in Tabled. 

This error bar only covers the random sources of error. To 
get a feeling for the systematic errors, we redid the reduction 
in a slightly different way. Instead of robust uniform weighting, 
we tried both uniform weighting (by setting the robust parame- 
ter to -2) and natural weighting {robust parameter = +2). In the 
latter, the RMS in the maps is lower at the expense of the more 
elongated beam and worse side-lobe levels. In the former, the 
beam shape is smaller at the expense of higher RMS levels in 
the map. The robust uniform and natural weighting flux deter- 
minations agree well within the error bar. As usual, the better 
the UV coverage and beam shape the less difference there is 
among the different flux determinations. We also measured the 
flux by determining the maximum intensity of the point source, 
instead of fitting an elliptical Gaussian. This alternative flux de- 
termination always falls within the error bar. 



3.4. Results 

The resulting 3.6 and 6 cm fluxes are plotted in the top panel of 
Fig. |3 While this figure suggests that £ Pup is variable at both 
wavelengths, we found a similar behaviour in other sources 
on the map. We therefore compared the flux values for £ Pup 
with two of the brightest sources close-by (S2=EQ 0801-398 
and S3=NVSS J080327-395828, see Fig. Q). The comparison 
(lower panels on Hig- EJ shows that all three sources follow the 
same trend. Hence, this argues strongly against £ Pup show- 
ing variability. A further argument against variability follows 
from the fluxes derived from the complete data set (see below) 
which are also plotted on Fig. [2] (dotted lines). Obviously, the 
flux of the complete data set is not the straight average of the 
separate fluxes. This points to problems in the interpolation of 
the instrumental gain phases during the first few runs of the ob- 
serving session, causing part of the £ Pup flux to be scattered 
over the rest of the image. These runs get less weight in the 
complete data set due to the non-linear nature of the cleaning 
process. Having concluded that the present data do not show 
detectable variability in £ Pup, we can derive upper limits on 
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the amount of variability from the range in fluxes (with their 
error bars): these are +20 % of its flux value at 3.6 and at 6 cm. 

In principle, the non-detection of variability could be 
caused by a bad coverage of the phase for the various periods. 
However, the 8 observations spread over 12 days cover the 5.2 
day rotation period reasonably well. We also checked the phase 
coverage of the variability cycle for possible periods around 8.5 
h and 19 h (see Sect. [2) and found it to be good. 

As there is no variability, we combined all visibility data 
and Fourier transformed them to obtain a single map. Because 
the RMS in this combined map is lower than in each individual 
map, we cleaned it deeper, down to 0. 1 mJy. The resulting im- 
ages at 3.6 and 6 cm are shown in Fig. [3] The measured fluxes 
and their error bars (RMS + 2% of the flux) are listed in Table|5] 

To get a feeling for the systematic errors on the combined 
data set, we investigated the effect of the cutoff limit in the 
CLEANing procedure: we CLEANed the total map at 6 cm (ro- 
bust-0.5, rms=0.024 mJy) down to 0.07 mJy and 0.15 mJy. 
In the former case the number of clean components was 2769 
in the latter 160. The resulting ( Pup fluxes fall well within the 
error bar. Also using a single, large, cleaning box instead of 
multiple cleaning boxes gave negligible differences in the flux. 
Finally, we redid the reduction, systematically dropping one of 
the observing runs. The range of values thus obtained again 
falls within the error bar. 



4. VLA observations 

4.1. 20 cm observation 

We obtained a 20 cm (L-band) observation on 2002 March 27 
with the VLA in A configuration (i.e. the configuration with the 
highest spatial resolution). The observation alternated between 
£ Pup and the phase calibrator 0814-356 (J2000). Two runs of 
7 min were made on ( Pup and 3 runs of 2 min on the phase 
calibrator. For the flux calibration we also observed 3C48 = 
0137+331 (J2000). The observation consists of two sidebands 
(at 1.3851 and 1.4649 GHz), each of 50 MHz bandwidth. 

The reduction of these data was done using the NRAO 
package AIPS (Astronomical Image Processing System), fol- 
lowing the same steps as for the ATCA data (Sect. [3}- Technical 
details of the reduction are listed in Table[5] We stopped clean- 
ing when the algorithm started finding about the same number 
of negative as positive components. The resulting map is shown 
in Fig. H 

By fitting an elliptical Gaussian to £ Pup, we found a 20 
cm flux of 0.76 ± 0.09 mJy. The error bar is the RMS noise 
in the total map. We also added a 2 % calibration error to the 
result (Perley & Taylor 2002 1, but this does not change the re- 
sulting error bar significantly. The maximum intensity of ^ Pup 
is 0.74 mJy/beam: the close agreement with the value from the 
Gaussian fit shows that £ Pup is indeed a point source. 

To judge the robustness of our flux determination, we re- 
peated the reduction, systematically dropping one antenna, 
using different weightings of distant visibilities, or natural 
weighting instead of robust uniform, or doubling or halving the 
number of clean components. In all cases, the results fall within 



AB1017 2002-MAR-27 20 cm 

1 ran rrn I i i nn 




b i i [ i i i i I i .llj i L 

08 h 03 m 36 s 35 s 34 s 

RA (J2000) 
-0.2,0.2,0.4,0.6 mJy 



Fig. 4. VLA observation at 20 cm. The contour levels follow a 
linear scale and the values are listed at the bottom of the fig- 
ure. The negative contour is given by the dashed line. The first 
positive contour is at about twice the RMS noise in the map. 
The beam is shown in the upper left corner. The cross indicates 
the optical ICRS 2000.0 position (from SIMBAD), corrected 
for proper motion. The ~ 1'.'4 offset with the radio position is 
most probably due to ionospheric refraction, which is more im- 
portant at longer wavelengths and at the low elevation at which 
f Pup was observed. 

the error bar. The RMS error bar therefore not only covers the 
statistical uncertainty, but the systematic errors as well. 

4.2. Archive observations 

The VLA archive contains a number of f Pup observations, 
which are listed in Tabled Many of these observations have 
not been published previously. To avoid introducing systematic 
effects we decided to reduce the whole set. The procedure is the 
same as that followed in Sect. 14. II The details of the reductions 
are summarized in Table [5] In principle, the 1.3 cm (K-band) 
and 2 cm (U-band) observations should be corrected for differ- 
ential atmospheric extinction (using the AIPS task elint), but 
in none of those cases was there sufficient information to apply 
this. 

The resulting maps are shown in Fig. [5] and the fluxes are 
listed in Table[5] The error bars show the range of results found 
by various reductions where we systematically dropped one an- 
tenna, used different weightings of distant visibilities, or nat- 
ural weighting instead of robust uniform. The error bars are 
always larger than the RMS noise in the map. For each obser- 
vation, we compared the peak intensity of the source to the flux 
derived from the Gaussian fit. With the only exception of the 
3.6 cm observation of AH365 (see below), there is very good 
agreement, showing that the sources are indeed point sources. 
When we doubled or halved the number of clean components, 
we were always well within the error bar. The error bars of the 
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Table 5. Reduction of the VLA data. Column (1) gives the wavelength, (2) the programme name, (3) the date of the observation, 
(4) the configuration the VLA was in at the time of the observation, (5) the name of the flux calibrator used, (6) the flux (in Jy) 
assigned to the flux calibrator (two values are given if there are two sidebands), (7) the phase calibrator (J2000 coordinates), (8) 
the number of antennas that gave a usable signal, (9) the integration time (in minutes) on the target, (10) the beamsize (major x 
minor axis) in arcsec 2 , (11) the measured flux (in mJy) and (12) refers to the notes. The fluxes of the flux calibrators are based 
on Perley & Taylor ( 20(1). 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


,1 


progr. 


date 


conf. 


flux calibrator 


phase 


no. 


intgr. 


beamsize 


flux 


notes 


(cm) 








name 


flux (Jy) 


calibrator 


ants. 


time 


(arcsec ) 


(mJy) 




1.3 


FLOR 


1981-10-18 


C 


3C286 


2.517 


0828-375 


19 


56 


5.6x 3.0 


< 13 


1,2 


2 


NEWE 


1979-02-09 




— 




0730-116 


8 


9 




< 13 


1,2,3 




FLOR 


1981-10-18 


C 


3C286 


3.452 


0828-375 


24 


64 


5.1x 3.5 


4.3 ± 0.9 


1 




AA28 


1984-03-07 


CnB 


JV.ZOD 


J.^+Z.J/J.HJZ. 


0828-375 


24 


18 


1.7x 1.3 


2.9 ± 0.3 


1.4 




AB327 


1985-01-29 


A 


3C48 


1.742/1.748 


0828-375 


26 


121 


3.0x 1.8 




1,5,6 


3.6 


AH365 


1989-05-13 


CnB 


3C48 


3.171/3.153 


0828-375 


26 


14 


4.0x 3.0 


1.4 ■ " i 


6.7 


6 


FLOR 


1978-07-23 




3C286 


7.462 


0836-202 


11 


9 




< 4 


8 




FLOR 


1978-10-13 




3C286 


7.462 


0836-202 


9 


160 


9.4x 1.3 


1.7 ±0.3 


9 




BIEG 


1978-11-05 




3C286 


7.462 


0836-202 


8 


24 


16. x 0.69 


< 4 


8 




CHUR 


1979-07-12 




3C286 


7.462 


0828-375 


12 


140 


4.9x 0.55 


1.4 ±0.3 


10,11 




FLOR 


1981-10-18 


C 


3C286 


7.462 


0828-375 


27 


28 


18. x 3.6 


1.71 ± 0.14 






AA28 


1984-03-07 


CnB 


3C286 


7.462/7.510 


0828-375 


24 


18 


5.0x 4.1 


1.49 ± 0.11 


12 




AB327 


1985-01-29 


A 


3C48 


5.405/5.459 


0828-375 


27 


91 


l.lx 0.45 


1.05 ±0.08 


6 


20 


FLOR 


1979-02-16 




3C286 


14.51 


0836-202 


11 


182 


14. x 2.3 


< 1.5 


2 




FLOR 


1981-10-18 


C 


3C286 


14.51 


0828-375 


26 


18 


60. xl2. 


< 0.75 


2 




AB1017 


2002-03-27 


A 


3C48 


15.49/16.20 


0814-356 


26 


14 


5.5x 1.1 


0.76 ± 0.09 





Notes: 

1 Insufficient data for atmospheric extinction correction 

2 Upper limit from 3 sigma 

3 No flux calibrator, used 0532+075 (J2000) instead with an assumed flux of 2.20 Jy 

4 Bieging et al. 1 1989 1 list 3.0±0.2 mJy for this observation 

5 Use of 3C48 as a flux calibrator in the A configuration is not recommended (Perley & Taylor 2002 1 

6 Data reduction problems detailed in Sect. 14. 21 

7 Lamers & Leitherer 1 1993 1 list 1.60±0.07 mJy for this observation, based on work by Howarth & Brown 119911 

8 Upper limit derived from non-detection of EQ 0801-398 

9 Primary is resolved; true flux might be slightly lower 

10 An additional 10 % uncertainty in the flux calibration should be added, because the instrumental gains of the flux and phase calibrator 
show significant differences (Fomalont & Perley 1999 1 

1 1 Abbott et al. 119801 and Bieging et al. 119891 list 1.4±0.3 mJy for this observation 

12 Bieging et al. 1 1989 1 list 1.3±0. 1 mJy for this observation 

3.6, 6 and 20 cm observations include a 2 % calibration error, 
and the 2 cm a 5 % calibration error. 

As ( Pup is rather low on the horizon from the VLA, the 
quality of the observations is not always good. The AB327 
observations at 2 and 6 cm show large variations in the gain 
phases (30-50 degrees) of the calibrators. All sources on the 6 
cm image (including ( Pup) are 30-50 % lower in flux com- 
pared to other observations. The 2 cm £ Pup flux is heavily de- 
pendent on tapering, and is compatible with a value of ~ 3 mJy. 
The 3.6 cm observation of AH365 is the only one where the 
measured peak intensity is systematically higher (by ~ 40 %) 
than the flux derived from the Gaussian fit. This suggests that, 
due to phase problems, the source is not quite a point source. 



The measured peak intensity is 2 mJy. Minor comments on 
other observations are given in Table[5] 

In general, the agreement of the fluxes we derived with pub- 
lished values is quite good (see notes to Tabled. Only for the 
AA28 - 6 cm observation do we find a value (1.49 + 0.1 mJy) 
significantly higher than the published value (1.3 + 0.1 mJy, 
Bieging et al. lT989l 

There is an additional observation available (AC308 - 
see Tabled made at 20 cm. This observation is part of the 
NRAO VLA Sky Survey (NVSS - Condon et al. 19981. We 
checked this survey and found that £ Pup was not detected (but 
EQ 0801-398 is detected). 
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Fig. 5. Same as Fig.0] but for the VLA archive data at 2, 3.6 and 6 cm. Only those observations having a detection are shown. 
Note the different spatial scales and different contour levels of the maps. 



4.3. Long term variability 

We now compare the different VLA observations, wavelength 
by wavelength, to see if we can detect variability at timescales 
much longer than the rotation period. 

We have only two 2 cm flux determinations. The error 
bars do not overlap, which might suggest variability. However, 
we recall the considerable difficulties in the reduction of 
these data (Sect. 14. 2> . and based on that, we do not consider 
them to present evidence of variability. For completeness, we 
also mention the 2 cm observations of £ Pup by Morton & 
Wright (119791 1: they list 7.2+1.1 mjy as the average of two 



observing runs. This value is significantly higher than ours. 
However, these observations were made with a single-dish an- 
tenna (Parkes 64-m), which has a beam of 2'3. When target- 
ing ( Pup, the beam also covers S2=EQ 0801-398 (see Fig.^, 
which will dominate the measured flux. These data can there- 
fore not be used to look for variability of £ Pup. 

The single VLA 3.6 cm observation (1.4 + 0.3 mjy) can 
be compared to the ATCA determination (2.38 + 0.09 mjy). 
In Sect. 14.21 we showed that the AH365 observation was 
quite problematic. Comparing the ATCA measurements of 
sources S2=EQ 0801-398 and S3=NVSS J080327-395828 to 
the AH365 measurements shows only a 20 % effect for S2, 
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Fig. 6. The 6 cm fluxes of £ Pup as a function of time. The open circle indicates an observation (AB327) that might be unreliable 
(see Sect. 14. 21 . The dotted line shows the flux from the ATCA combined data. 



and none for S3, suggesting that the Gaussian fit fluxes would 
be reliable. In that case, the AH365 £ Pup flux would be sig- 
nificantly different from the ATCA determination. However 
AH365 shows such an accumulation of problems (the flux of 
the secondary calibrator is considerably less than expected, 
there are reasonably high phase changes and there is the dif- 
ference between peak intensity and Gaussian integrated flux) 
that we do not consider this convincing evidence of variability. 

The single 20 cm determination is compatible with the two 
upper limits. 

The only wavelength for which we have a reasonable num- 
ber of flux determinations, is 6 cm. Fig. [6] shows these fluxes 
as a function of time. The flux from the ATCA combined data 
goes through the error bars of almost all observations, indicat- 
ing that there is no detectable variability. 

The only exception is the significantly lower flux of AB327. 
In Sect. 14.21 we discussed the phase problems of this obser- 
vation, which lead to the lower flux. Another possible factor 
might be that we are starting to resolve the stellar wind of 
( Pup (this is the observation with the highest spatial resolu- 
tion we have). To estimate how much flux we would lose, we 
use a Wright & Barlow (119751 approach. By their definition of 
the characteristic radius (R v ) of the radio-emitting region, the 
observed flux (S v ) is given by: 

v 2 



S v = 10 



-26 



f 



dr-^K{v, T)y 



M 



47IT 2 l 



B V {T), (1) 



where r is the radius, D the distance to the star, K(v, T) the free- 
free absorption coefficient, y the ratio of electron to ion number 
densities, M the mass loss rate, Voo the terminal velocity, fi the 
average atomic mass, mu the proton mass and B V (T) the Planck 
function at frequency v and temperature T. All units are cgs, 
except for S v which is in mJy. If, instead of integrating to +oo, 
we only integrate to a radius R max , we lose a fraction of the flux 
given by: 

S v ^max 



The beamsize (#b eain ) determines fl max by: 6» be am = 2R max /D. 
Combining this with the Wright & Barlow expression for R v , 
we get: 



AS, 
S v 



3.8 x 10 : 



XygZ 2 y /3 T 



-1/2 



MA 

jUVoo 



2/3 



(3) 



where M is in Mg/yr, v M in km/s, A in cm, #b eam in arcsec, D 
in kpc. We take M = 5.9 x 10~ 6 (Puis et al. fT99%t . fi=lA, y=l, 
Z=l and the Gaunt-factor g-1 . Other parameters are taken 
from Table 13 With a beam of 6*b e am=077, we find a flux loss 
of 6 - 12 %, depending on whether we take a hot wind (42 000 
K) or a cool wind (10000 K). If we use the mass loss rate de- 
rived from our model for the radio observations (Table 0, we 
find a 4 - 8 % effect. 

The low AB327 flux can thus not be explained by our start- 
ing to resolve the stellar wind, but is due to problems with the 
interpolation of the gain phases. 

5. JCMT observation 

We also determined the 850 yum flux of £ Pup using 
the Submillimetre Common-User Bolometer Array (SCUBA, 
Holland et al. [T9^9l on JCMT (see Table [I). The instrument 
observes simultaneously at 450 and 850 fim using two hexag- 
onal arrays of bolometers. The sensitivity at 450 fim is too low 
for a detection, so we will only discuss the 850 yum data. Our 
flexibly scheduled observations were taken in medium weather 
conditions (850 fim zenith opacity was 0.32 - 0.47) during a 
half-shift on 2000 October 12. The beam at 850 fxm is 1475. 

As the data were collected in the same run as the e Ori 
observation discussed by Blomme et al. (2002i, we refer to 
that paper for details on the observation and the reduction. The 
total on-target integration time for £ Pup was 15 min. After 
reduction, the RMS on the £ Pup observation turns out to be 
about 16 %. We tried variants in the reduction (e.g. using other 
bolometers than the inner ring for sky-noise removal), but this 
changes the flux by considerably less than 16 %. 
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Table 6. Submillimetre and radio fluxes of ( Pup. Data from 
this paper, unless otherwise indicated. 



A 


flux (mjy) 


reference 


850 /an 


28 ±5 


JCMT 


1.3 mm 


20.2 ± 1.8 


from Leitherer & Robert I1991i 


2 cm 


4.3 + 0.9 


VLA, FLOR 1981-10-18 


2 cm 


2.9 ± 0.3 


Bieging et al. ( 1989 1, revised 


3.6 cm 


2.38 ± 0.09 


ATCA 


6 cm 


1.64 ±0.07 


ATCA 


20 cm 


0.76 ± 0.09 


VLA 



The £ Pup observations are preceded and succeeded by an 
observation of the calibrator OH231.8. Taking the average of 
these two calibration observations, we arrive at 31 + 5 mjy for 
the flux. The error bar takes into account the measurement er- 
rors on the target and the calibrator as well as the calibration 
error. However, OH231.8 is somewhat variable and the flux of 
the calibrator at the time of our observation could be different 
from the reference flux we used. 

If we use the HL Tau calibration observation that precedes 
our run, we have 28+5 mjy. Trying to use other calibrators that 
were observed during that night (but of course further away in 
time from our observation) tends to favour the lower value. In 
view of this we propose 28+5 mjy as the best determination of 
the flux. The ( Pup value is then also based on the same calibra- 
tor as the e Ori observation analysed by Blomme et al. (2002 1. 

If we want to compare our 850 /mi value with the 1.3 mm 
determination of Leitherer & Robert d!991> . we need to cor- 
rect for the difference in wavelength. If we assume an a — 0.6 
spectrum, we find that our value corresponds to 22 + 4 mjy at 
1.3 mm, which agrees very well with the Leitherer & Robert 
value of 20.2 ±1.8 mjy. Our value is also compatible with the 
Altenhoff et al. ( 1994 1 upper limit of 33 mjy at 1.2 mm. 

6. Discussion 

6.1. Smooth wind model 

Table[6]summarizes the submillimetre and radio fluxes we will 
discuss in this section. Other fluxes either agree with these, or 
possible disagreements have been explained in Sect. I4.2I 

From the 1.3 mm and 6 cm fluxes we can derive the spec- 
tral index a (defined by F v oc A~ a ). The measured value of 
a = 0.66 + 0.03 defines a power law that goes through the error 
bars of all observations, with the exception of the 2 cm fluxes 
(where it is intermediate between the two determinations). This 
spectrum is somewhat steeper than expected from the Wright 
& Barlow ( 1975 1 model (where a = 0.6), indicating some dis- 
crepancy between the millimetre and radio fluxes. 

To better quantify this discrepancy, we also made a smooth 
wind model for ( Pup in the same way as in Runacres & 



Blomme (1996 1. The model solves the equations of radiative 
transfer and statistical equilibrium in a spherically symmetric 
stellar wind, containing only hydrogen and helium. The den- 
sity in the wind is determined by solving the time-independent 
hydrodynamical equations (following Pauldrach et al. [1986 1. 
When fitting the model to the observations, the visual and 




lOOnm 1 fj,m 



10 (MTn 100 /um 1 mm 
wavelength 



1 cm 10 cm 



Fig. 7. Observed fluxes normalised to a smooth wind model. 
Observations above the dotted line point to additional emission 
that is not included in the smooth model. 

Table 7. Comparison to other mass loss rate determinations of 
<TPup. 



wavelength M(10 MqIjx) 
9 4+1-0 

3.5 +L9 



reference 



radio 
Ha 



radio 



5.9 
3.5 



-1.2 



Lamers & Leitherer 1 1993 1 
Lamers & Leitherer 1 1993 1 
Puis et al. fT996l 
this paper 



near-infrared fluxes were used to determine the interstellar ex- 
tinction and the radio fluxes to determine the mass loss rate. 
The far-infrared and millimetre fluxes are unconstrained and 
can therefore be used to see how well the smooth wind model 
fits the observations. Further details of the model are given in 
Runacres & Blomme. A new version of their £ Pup model was 
calculated that uses the parameters listed in Table[2] 

In fitting the model, we fixed the mass loss rate from our 
ATCA 6 cm observation. We choose this observation because 
it has the smallest error bar. The best fit gives a mass loss rate of 
M = 3.5 ± 0.1 x 10~ 6 M Q /yr. The error bar on M corresponds 
only to the error bar on the 6 cm flux, and does not include 
the more important errors due to stellar parameters or distance. 
E.g., the error in the Hipparcos distance converts to a (-0.9, 
+ 1.5) xlO" 6 M Q /yr error on M. 

Because it is such a well-studied star, a large number of 
mass loss rate determinations have been made for £ Pup. In 
our comparisons we will only consider the more recent deter- 
minations (Table 0. Our radio mass loss rate of 3.5 x 10 -6 
Mo/yr is different from the 5.9 x 10 -6 M©/yr value by Puis 
et al. ( 1996 1 based on fitting the Ha line profile. This differ- 
ence is significant, considering that we used the same stellar 



parameters as they did. We recall that Petrenz & Puis (1996 1 
already noted this discrepancy and tried to explain it by invok- 
ing an equatorial density enhancement, caused by the large ro- 
tational velocity of the star. Our results are in acceptable agree- 
ment with both the radio and Ha determinations by Lamers & 
Leitherer ( 1993 1, but we note that they used a simpler model for 
the Ho- emission than Puis et al. Contrary to what the Lamers 
& Leitherer numbers suggest, £ Pup is therefore an exception 
to the usually good agreement between the Ha and radio mass 
loss rates. 

The observed fluxes are compared to the best-fit smooth 
wind model in Fig. [7] The figure shows the millimetre and ra- 
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Fig. 8. The millimetre excess is fitted with a structured wind model. Three possible runs of the clumping factor as a function of 
radius are used (inset). The dashed line presents a model that has too much structure too far out in the wind to be able to explain 
the observations. 



dio observations listed in Table0 supplemented by visual and 
infrared data listed in Runacres & Blomme ( 1996 1. As found 
by Runacres & Blomme, there is an excess flux at millimetre 
wavelengths. Due to our somewhat higher 6 cm flux, the effect 
found here (24 %) is somewhat lower than what they found 
(35 %). The present result is significant at the 2 sigma level. In 
view of the large error bars, it is not clear what happens at 2 cm, 
but the 3.6 cm flux could be marginally in excess. The 20 cm 
flux is in good agreement with the smooth wind model. We 
stress that the same discrepancy is found when a much simpler 
Wright & Barlow ( i!975i model is used instead of our smooth 
wind model. 

6.2. Explaining the millimetre excess 

The millimetre excess is the most important result from this 
work. Various explanations will be considered in this section. 
To study the excess we will make considerable use of models 
similar to the one developed by Wright & Barlow jl975> . 

We start by noting that the Wright & Barlow formula for 
the flux is not very sensitive to temperature, so a radial gra- 
dient of the temperature will not explain the excess. There 
could be an indirect effect of the temperature however, due 
to the recombination of important ions. This explanation has 
been proposed by Leitherer & Robert ( 1199U . As He ++ will re- 
combine much more quickly than H + , we only consider the 
recombination of He ++ . A careful evaluation of the ionization- 
dependent factors in the Wright & Barlow dl975> formula for 
the radio flux shows that, if we assume all helium to be He ++ 
in the millimetre formation region and all helium to be He + 
in the centimetre formation region, the discrepancy between 
the fluxes can be explained. It should be noted that our own 
smooth-wind model does not show recombination, but it does 
not include important effects, such as line blanketing due to 



metals. Recent models of early-type stars arrive at lower ef- 
fective temperatures than the one we used here (Bianchi & 
Garcia 120021 Bianchi et al.[2003 Martins et al. l2U02l . which 
would favour recombination. Further evidence that recombina- 
tion occurs, comes from models of the stellar wind of ( Pup 
(Hillier et al. [Tff^l Pauldrach et al. I2TKJT1 Puis 2003, pers. 
comm.). Due to the recombination in these models, it is neces- 
sary to assume that the X-rays are formed out to large distances 
in the wind (> 100 R t ), because X-rays formed too close to the 
star get absorbed. However, the X-ray spectral lines of f Pup 
(Cassinelli et al. 2001 Kahn et al. 2001 1 seem to require for- 
mation closer to the stellar surface and less opacity in the wind 
(Kramer et al. 2003 1. This fact is difficult to reconcile with He 
recombination. 

A slower velocity law has been claimed to explain the far- 
infrared discrepancies seen on Fig.0(Kudritzki & Puis 2000 1. 
To explain the millimetre discrepancies would require a very 
much slower velocity law. To investigate this we adapted the 
Wright & Barlow (1975) model to include a ft-type velocity 
law. In such case, an analytical solution is no longer possible, 
and all integrations were done numerically. We find that ft ~ 5 
is needed to explain the millimetre fluxes, which is not compat- 
ible with other observational indicators such as the Ha profile, 
which requires ft ~ 1.15 (Puis et al. 1996 1. 

The existence of a flux excess at millimetre wavelengths 
can also be due to structure in the wind. As there is a wealth 
of evidence for structure in these winds, it seems obvious to 
explore a model based on that. 

The model in question is the same as used in Blomme et 
al. (|2002) for e Ori. The smooth wind opacity in a Wright & 
Barlow l !1975t model is multiplied by the clumping factor to 
get the clumped wind opacity. The clumping factor is defined as 
< p 2 >/< p > 2 , where the symbol <> stands for a time-average, 
which we have approximated by integrating over a small vol- 
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ume of the wind. Optical depth, emergent intensity and flux are 
then calculated as in Wright & Barlow. As our model allows 
the clumping factor and velocity to change as a function of dis- 
tance, all integrations have to be done numerically. We used a 
run of the clumping factor based on the work of Runacres & 
Owocki (2001 1, who calculated time-dependent hydrodynam- 
ical models to study the effect of the line-driving instability 
at large distances from the star. The clumping factor in their 
models rises to reach a maximum rather far away from the star 
(10 - 50 R*) and then decreases again. If we simplify their re- 
sults somewhat, we can approximate the typical behaviour for 
the clumping factor by a piece-wise linear curve fixed by speci- 
fying three points. We let the clumping factor rise linearly from 
close to the surface of the star (point 1) to a certain distance 
(point 2), and then let it fall off again linearly till it becomes one 
(point 3). The inset to Fig.[8]shows an example (dotted line). 
We also used piece-wise linear curves fixed by four points. By 
setting the clumping factor equal to one at large distances, we 
ensure that there is no excess at radio wavelengths. 

Fig- IH shows a number of fits we attempted. Although no 
unique solution is found, it is quite clear that the clumping fac- 
tor has to be significantly higher than 1 in the inner part of the 
wind, and needs to diminish considerably beyond 70 R*. To 
stress this last point, Fig. [8] includes a model with substantial 
clumping up to 120 R* (dashed line); the fluxes clearly over- 
shoot the 3.6 and 6 cm observations. The fact that structure di- 
minishes beyond a certain distance is similar to what we found 
for e Ori (Blomme et al. 2002 1, except that for e Ori, this started 
happening around ~ 40 R,. 

While structure diminishes around 70 R t , it does not nec- 
essarily disappear. For e Ori we showed that a model where 
structure persists up to large distances (with a constant clump- 
ing factor) could equally well explain the observations, pro- 
vided we reduce the mass loss rate accordingly. For f Pup, we 
see that the 20 cm observation is in reasonable agreement with 
the 6 cm one. This shows that, there is still clumping left, 
it falls off considerably slower than it does in the inner 70 R* 
of the wind. Obviously, a more accurate 20 cm flux and obser- 
vations in the 1 mm - 3.6 cm region will further constrain the 
extent and amount of structure. 

The above model was very much inspired by the small- 
scale structure expected from the instability of the radiative 
driving mechanism. However, CIRs, a disk or a polar enhance- 
ment should have similar effects on the millimetre flux. In these 
cases the 70 R, radius will be significant as well, in that it 
shows where the structure diminishes substantially, or maybe 
even disappears. A detailed comparison of models for different 
types of structure is beyond the scope of the present paper. 

7. Conclusions 

Radio observations of f Pup covering about two rotational pe- 
riods, supplemented by archive observations covering a much 
longer time scale, do not show variability at more than the 
+20 % level. The long integration time gives us an accurate 
flux determination of 2.38 + 0.09 mJy at 3.6 cm and 1.64 + 
0.07 mJy at 6 cm. These values are slightly higher than the 
ones previously known. 



Converting the fluxes into a mass loss rate, we find M = 
3.5 x 10~ 6 Mq/jt. This value confirms the significant discrep- 
ancy with the Ha mass loss rate (Petrenz & Puis 1996 1. 

A smooth wind model shows that the millimetre fluxes are 
too high compared to the radio fluxes. While recombination 
of helium in the outer wind cannot be discounted as an ex- 
planation, we favour a model that ascribes the discrepancy to 
structure. A simple model shows a substantial decay, or maybe 
disappearance, of structure beyond 70 R,. Fig.[8]shows how ob- 
servations at wavelengths between 1 mm and 3.6 cm can further 
constrain models for structure. 

The present data do not allow a distinction between the 
various types of structure (stochastic, CIRs, disk or polar en- 
hancement). Attempting to detect variability at far-infrared or 
millimetre wavelengths should provide better constraints on the 
azimuthal symmetry of the structure, thereby allowing us to de- 
cide whether it is in the form of CIRs. 
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